Solving the Graph Isomorphism Problem with a Quantum Annealer 
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We propose a novel method using a quantum annealer - an analog quantum computer based on 
the principles of quantum adiabatic evolution - to solve the Graph Isomorphism problem, in which 
one has to determine whether two graphs are isomorphic (i.e., can be transformed into each other 
simply by a relabeling of the vertices). We demonstrate the capabilities of the method by analyzing 
several types of graph families, focusing on graphs with particularly high symmetry called strongly 
regular graphs (SRG's). We also show that our method is applicable, within certain limitations, to 
currently available quantum hardware such as "D-Wave One". 



I. INTRODUCTION 



Theoretical research on quantum computing is motivated 
by the exciting possibility that quantum computers may 
be able to perform certain tasks faster than classical 
computers. In recent years first steps have been taken 
towards the goal of experimentally realizing these com- 
putational advantages. However, to date the largest 
experimental implementations of scientifically meaning- 
ful quantum algorithms have used just a handful of 
qubits 0,0, the reason being the tremendous technolog- 
ical challenges (the most crucial of which is overcoming 
quantum decoherence) that need to be defeated before 
a successful implementation of any solid-state quantum 
computer. 

Of the several quantum computing paradigms that have 
been proposed, a potentially promising substitute for the 
'standard' circuit-based quantum computer is the quan- 
tum annealer 0, Q , which is now on the cusp [|| of be- 
ing able to run small-scale computing procedures on ac- 
tual quantum annealing hardware. Quantum annealing 
machines are analog quantum computational devices de- 
signed to solve discrete combinatorial optimization prob- 
lems using properties of quantum adiabatic evolution. 
They are based on a general approach widely known as 
the Quantum Adiabatic Algorithm (QAA), which was 
proposed by Farhi et al. [5| about a decade ago as a 
method for solving a broad range of optimization prob- 
lems using a quantum computer. 

Recently, a quantum annealing machine, based on liq- 
uid crystal nuclear magnetic resonance has been reported 
to successfully factor the number 143 using four spin- 
qubits @. A few months later, "D-Wave One", a 128- 
qubit machine based on super-conducting qubits Q ex- 
perimentally demonstrated the ability to compute two- 
color Ramsey numbers @, H[ making it the largest ex- 
perimental implementation of a scientifically meaningful 
quantum algorithm. This success provides an incentive 
for finding problems that could be solved efficiently on a 
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quantum annealer. 

In this paper, we hypothesize that a quantum annealer 
could solve the Graph Isomorphism (GI) problem, de- 
scribed in detail below. In particular, we hypothesize 
that the annealer could distinguish all non-isomorphic 
graphs by sufficiently precise measurements. We pro- 
vide some evidence for this by showing that our method 
works for certain graphs which we have been able to 
study numerically. An advantage of our proposal is that 
it can, within certain limitations, be implemented using 
currently available quantum annealing machines such as 
D-Wave One. 

The paper is organized as follows. In Sec.|TT]we describe 
the Graph Isomorphism (GI) problem in some detail. 
In Sec. |nT] we describe and then discuss the proposed 
method for solving the GI problem using quantum adia- 
batic evolution. We next present the results of the study 
in Sec. II VI and derive some conclusions in Sec. fVl 



II. THE GRAPH ISOMORPHISM PROBLEM 



A graph G = (V, E) is a set of vertices V, and edges E 
which arc unordered pairs of vertices. A graph is con- 
veniently expressed algebraically as an adjacency matrix. 
The adjacency matrix A of a graph with N vertices is an 
N x N matrix in the basis of vertex labels, with Ay = 1 
if vertices i and j are connected by an edge, and zero 
otherwise. 

The GI problem is stated as follows: Given two graphs, 
one must determine whether or not they arc isomorphic 
to each other, i.e., whether one can be transformed into 
the other by a relabeling of the vertices. 



In the realm of classical computing, many special cases of 
GI have been shown to be solvable in a time that scales as 
a polynomial of the number of vertices. However the best 

general algorithm to date runs in time O (^a N ' lo & N \ 

where a is a constant Q . It is therefore interesting to ask 
whether a quantum computer could solve this problem 
efficiently. An attractive feature of the GI problem is 
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its resemblance to the problem of integer factoring, the 
first and most famous example to date of a quantum 
algorithm can solve a problem exponentially faster than 
the best known classical algorithm Like factoring, 

the common belief is that the GI problem is unlikely to 
be NP-complete. 

The GI problem has been attacked by numerous meth- 
ods inspired by classical as well as quantum physical sys- 
tems. Rudolph mapped the GI problem onto a sys- 
tem of hard-core atoms, where one atom was used per 
vertex, and atoms i and j interacted if vertices i and 
j were connected by edges. It was demonstrated that 
for some pairs of non-isomorphic graphs, sharing cospec- 
tral adjacency matrices does not lead to cospcctrality of 
the transition matrix between three-particle states pro- 
ducedby the embedded Hamiltonian. Gudkov and Nussi- 
nov [12( proposed a classical algorithm to distinguish 
non-isomorphic graphs by mapping them onto various 
physical problems. Shiau et al. proved that the simplest 
classical algorithm fails to distinguish some pairs of non- 
isomorphic graphs and also proved that continuous-time 
one-particle quantum random walks cannot distinguish 
some non-isomorphic graphs [l3l - [T5j . More recently, it 
has been found that classical random walks and quantum 
random walks can exhibit qualitatively different proper- 
ties (l6l - fl8l |. These disparities mean that in some cases 
algorithms implemented by quantum random walkers can 
be proven to run faster than the fastest possible classical 
algorithm [THH. 



III. QUANTUM ADIABATIC EVOLUTION 
AND THE GI PROBLEM 



To solve the GI problem using quantum adiabatic evolu- 
tion, we assign to each graph G the Hamiltonian H{G), 
where 

H{G) = (1 - s)H d + sH p (G) , (1) 

which depends on a parameter s. For s = 0, H(G) is 
the standard 'driver' Hamiltonian for QAA algorithms, 
namely 

1 N 

H d = -J2°?, (2) 

i 

i.e., a transverse- field Hamiltonian, while for s = 1, H{G) 
is the 'problem' Hamiltonian H P (G) which is constructed 
according to the topology of the graph. A simple plausi- 
ble choice for H p (G) is 

H P (G) = ]T ofo?, (3) 
(«>eG 

i.e., an Ising antifcrromagnet on the edges of the graph. 
While the driver Hamiltonian is diagonal in the Y[ °~ x 



basis, the problem Hamiltonian is diagonal in the Y[ °~ z 
basis. 

The system is first prepared in the ground state of the 
driver Hamiltonian Hd, which is straightforward. The 
adiabatic parameter s is then varied slowly and smoothly 
with time from to 1, so that the Hamiltonian is con- 
tinuously modified from Hd to H p (G). If this process is 
done slowly enough, the adiabatic theorem of Quantum 
Mechanics (see, e.g., Refs. [29[ and Ho[) ensures that the 
system will stay close to the ground state of the instan- 
taneous Hamiltonian throughout the evolution. 

The premise of the method we suggest here is that 
the state of the system along the adiabatic evolution, 
i.e., the instantaneous ground state of the Hamiltonian, 
Eq. (fTJ , stores enough information to reflect the complex 
structure of the graph-dependent problem Hamiltonian 
H p (G), and that carefully chosen measurements along 
the adiabatic path will provide enough information to 
differentiate non-isomorphic graphs. 

The choice of the problem Hamiltonian in Eq. ([3]) is, 
of course, only one possibility. However it has the ad- 
vantages of (i) using Ising spins which are simple to 
study and to implement experimentally @, Q, and (ii) 
having antiferromagnetic interactions which, on highly- 
connected graphs, tend to have highly frustrated ground- 
states because closed paths (along the edges) of odd 
length make the system a spin glass j3ll - l33 |. Note also, 
that the Hamiltonian Eq. (fTJ is symmetric with respect 
to flipping all spins. 

It is interesting to compare our approach with those 
which use Hamiltonians of quantum random walkers 
embedded in the graph structure. In the latter case 
the walker normally accesses only low-dimensional sub- 
spaces of the Hamiltonian eigenstates, based on con- 
served quantities such as number of particles (see, e.g., 
Refs. [ll|, [TH, H3, HH), whereas in our approach there is 
no such conserved quantity so the instantaneous ground 
state of the Hamiltonian presumably reflects the full com- 
plexity of the graph. 

For each graph, one performs multiple runs along the adi- 
abatic path, performing various measurements at differ- 
ent values of s until sufficient statistics is gathered. Since 
each measurement collapses the state of the system, non- 
commuting measurements or measurements correspond- 
ing to different s values, require separate runs of the pro- 
cedure. Since errors of the various measured quantities 
are inversely proportional to square root of the number 
of measurements, this number will be determined from 
the needed resolution. 

In order to make sure that isomorphic graphs are recog- 
nized as such, the measurements must be invariant un- 
der a relabeling of the vertices of the graph (or equiva- 
lently the spins in the system) . The most straightforward 
such quantity is the Hamiltonian, so the average energy 
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E = (H(G)) is a 'good' quantity to measure. Here, (■) 
indicates the expectation value with respect the state of 
the system, i.e., the instantaneous ground-state, at a par- 
ticular value of s. For the same reason, the classical 'di- 
agonal' average energy of the graph Eq = (H p (G)} and 
M x = 2(Hd), the x-magnetization, are also suitable ob- 
servables. 

Many other quantities that respect the topology of the 
graph and are invariant under any relabeling of the ver- 
tices could be measured. Some may prove to have better 
distinguishing capabilities than others, depending on the 
manner in which they 'tap' into the complexity of the 
structure of the graph in question. 

We find that a 'good' quantity to measure is the spin- 
glass order parameter, which we shall denote as Qi and 
define by 

Q2 =(]v(^g^ 2 ) • < 4 > 

Generalizing the above expression, the quantities 

Q^ = j^ E K<---<J 2 ; ( 5 ) 

where n — 2, 3, 4 . . ., also serve as distinguishing opera- 
tors. Analogously, other types of operators that might 
work are of the form: 

Q'2 = j!pIl(H p (G)a?cT*) 2 . (6) 

Each such measurable quantity, if measured sufficiently 
accurately, for different values of the adiabatic parame- 
ter s, serves as an additional 'dimension' of differentia- 
bility of the non-isomorphic graphs. Moreover, gather- 
ing statistics of these quantities for different values of s 
accesses indirectly the entire spectrum of H p (G) (unlike 
the situation in classical or quantum random walks where 
the walkers are restricted to only a small subspace of the 
Hamiltonian eigenstates) . 

In Sec. HVI we will consider several sets of non-isomorphic 
graphs, showing that, in all cases, non-isomorphic graphs 
are distinguished by accurate measurements of carefully 
chosen observablcs along the adiabatic path. In fact, we 
shall see that in most cases studying the limit s — > 1 
suffices. 



IV. RESULTS 



We now present numerical and semi-analytical results for 
several types of graph families that are known to be hard 
to distinguish. Since the annealing process requires that 



the temperature of the system be well below the excita- 
tion gap of the system, we shall work at zero temperature. 

For the smaller graphs that we study (N < 25 vertices), 
we use exact-diagonalization or conjugate-gradient based 
minimization techniques. Both techniques are restricted 
to quite small values of N because the size of the Hilbcrt 
space grows as 2 N . The specific variation of conjugate- 
gradient method that we developed for this study is dis- 
cussed in Appendix [X] 

For graphs of more than 25 vertices, quantum Monte- 
Carlo techniques would usually be the method of choice 
for the accurate measurements of the various quantities. 
However, these were found by us to be rather inefficient in 
the interesting regions where the value of the adiabatic 
parameter is close to 1. We have been able to study 
sizes a little larger than 25 (up to 29) by leading-order 
degenerate perturbation theory about the limit s = 1. 
This calculation involves finding the subspace of ground 
states of H p (G), which requires one to first calculate the 
energies of the 2^ 'classical' states followed by further 
manipulation on the subspace spanned by the ground- 
state eigenvectors. The main shortcoming of the method 
is that it only produces results in the s — > 1 limit. For- 
tunately, we find that in most cases this is sufficient. 

The measured quantities we focus on are the total (av- 
erage) energy, the spin-glass order parameter Q2 given 
in Eq. (Q| and the cc-magnetization M x = 2(Hd). We 
find that in most cases these three quantities are suf- 
ficient to distinguish all tested non-isomorphic graphs, 
though these measurements may be augmented by mea- 
suring other observablcs, as discussed in the previous sec- 
tion. 



A. Strongly regular graphs 

The main results of this paper are for strongly regular 
graphs (SRGs), a class of graphs, subsets of which are 
known to be difficult to distinguish [l4[. An SRG is a 
graph in which (i) all vertices have the same degree, (ii) 
each pair of neighboring vertices has the same number of 
shared neighbors, and (iii) each pair of non-neighboring 
vertices has the same number of shared neighbors. This 
definition permits SRGs to be categorized into families 
by four integers (N, k, A, fj,), each of which might contain 
many non-isomorphic members. Here, N is the number 
of vertices in each graph, k is the degree of each vertex (k- 
rcgularity) , A is the number of common neighbors shared 
by each pair of adjacent vertices, and \x is the number of 
common neighbors shared by each pair of non-adjacent 
vertices. 

Using the stringent constraints placed on SRGs, one can 
show that, for any SRG, the spectrum of the adjacency 
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FIG. 1: Spin-glass order parameter Q2 in the ground state 
for the two non-isomorphic strongly-regular graphs (SRGs) 
on TV = 16 vertices. As the figure indicates, the two graphs 
show different Q2 values starting from s ~ 0.4. 



matrix has only three distinct values [35|: 
which is non-degenerate, and Ai,2 = — \ ~ 



A = -k, 

which are both highly degenerate. Both the value and de- 
generacy of these eigenvalues depend only on the family 
parameters, so within a particular SRG family, all graphs 
are cospectral UK- These highly degenerate spectra are 
one reason why distinguishing non-isomorphic SRGs is 
difficult. 

While there exist SRG families with only one non- 
isomorphic member (36j , we concern ourselves with fam- 
ilies that have more than one non-isomorphic graph. Us- 
ing combinatorial techniques (37l . [38j , tables of complete 
and partial families of SRGs have been tabulated and we 
will use those tables to select graphs for our study. It 
should be noted that for every family of non-isomorphic 
SRGs there exists a complementary family of the same 
size N and the same number of members, which is ob- 
tained by interchanging edges and non-edges. 

The smallest family of non-isomorphic SRGs that have 
more than one member is that of the N = 16 vertices 
with signature (16,6,2,2), which contains two members 
(we shall not address here the complementary families of 
graphs which are distinguished in much the same way). 
The two graphs are immediately distinguished by looking 
for example at the spin-glass order parameter Q2 as a 
function of the adiabatic parameter s. The Q2 values 
of the two graphs (at zero temperature) are plotted as 
a function of s in Fig. [TJ In this N = 16 case, both 
graphs have the same 'classical' ground-state energy of 
Eq = (H p (G)) = —16 albeit with different degeneracies, 
namely 21 and 45 

The next family of SRGs that contains more than one 
member is that with N = 25 vertices. This family has 
15 distinct graphs with signature (25,12,5,6). We find 
that looking at the values of Q2 and M x in the ground 
state in the limit of s — > 1 using first-order degenerate 
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FIG. 2: Scatterplot in the M x — Q2 plane of the 15 strongly 
regular graphs with N = 25 vertices in the limit s — » 1. The 
horizontal axis is the magnetization along the i-direction and 
the vertical axis is the spin-glass order parameter Qi- The 
circled data point corresponds to the only two graphs which 
are not distinguished by values of Q2 and M x in the s — > 
1 limit. These two instances are, nonetheless, distinguished 
by measurements at s < 1 (see text). The inset shows the 
M x — Q2 values of two of the graphs that lie outside of the 
region shown in the main panel. 



perturbation theory suffices to distinguish between all 
but two of the graphs (the latter two share the same Q2 
and M x values in this limit). The results are shown in 
Fig. [21 which is a scatterplot of Q2 vs M x for ,s — > 1 for all 
15 non-isomorphic graphs. In this limit, most instances 
have the same ground-state energy of Eq = — 34 except 
for two instances which have a slightly higher energy of 
Eg = — 30. These correspond to the two points in the 
inset of Fig. [2] In addition, the circled data point in 
the figure corresponds to the (M X ,Q2) value of the two 
graphs which are not distinguished in the s 1 limit. 
However, for values of s less than 1, exact diagonalization 
reveals a clear distinction between the two graphs. For 
example, for s = 0.73 the values of Q2 in the ground 
state are 0.57914 and 0.443423. 

Perturbation theory-based analysis in the s — > 1 limit 
also allows us to study families of graphs with N = 26, 28 
and 29 vertices. The degenerate perturbation-theory 
analysis of the set of 10 SRGs with N = 26 and sig- 
nature (26,10,3,4) are shown in Fig. [3] which is a scat- 
terplot in the M x — Q2 plane. Here, 7 of the instances 
were found to have Eg = — 34 (with different degenera- 
cies, all of them around ~ 1000) and required first-order 
perturbation theory, whereas the remaining 3 instances 
had Eg — —38 (with degeneracies 20, 20 and 60), and the 
leading-order perturbation analysis is of degree 6. Fig. [3] 
shows that all 10 members of the family are distinguished 
in that plane. The inset shows data points that lie out- 
side of the range presented in the main panel. 

For N = 28, the family (28,12,6,4) contains 4 distinct 
graphs. There too, we find that it is sufficient to look at 
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FIG. 3: Scatterplot in the M x - Q 2 plane of the 10 
strongly regular graphs with N — 26 vertices and signature 
(26, 10,3,4), in the s — > 1 limit. As the figure indicates, the 
annealer distinguishes all graphs in the family. 
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FIG. 4: Scatterplot in the M x — Q2 plane of the 41 strongly 
regular graphs with N = 29 vertices, in the s — > 1 limit. As 
the figure indicates, the annealer distinguishes all graphs in 
the family. The inset shows data points that are out of the 
range presented in the main panel. 



the Q2 values and average total energy in the s — > 1 limit. 
For three of the four graphs that have the ground-state 
energy of Eq = —28 (and with different degeneracies, 
all of them around ~ 6000) the Q2 values are 0.137461, 
0.141957 and 0.132883. The fourth graph has a different 
energy of Eq = —24 (with a degeneracy of 972265). The 
Q2 value has not been calculated for this graph. 

The largest size of graphs that we deal with analytically 
using perturbation theory is the family of 41 graphs hav- 
ing N = 29 vertices and signature (29, 14,6, 7). A scat- 
terplot in the M x — Q2 plane for this family is given 
in Fig. 2J While some of the values shown in the scat- 
terplot lie close to one another, all pairs (M X ,Q2) are 
distinct (and can be further differentiated by measure- 
ments of other observables) so all non-isomorphic graphs 
of this family can be distinguished. Here, all members 
were found to have the same energy of Eq = —41 and 
with different degeneracies, all of them around ~ 3000. 
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FIG. 5: The difference in Q2 values, AQ2 (main panel) and 
energies (inset), for the two graphs Gi and G2 with N = 14 
that certain quantum walkers cannot distinguish [27I. l2fij]. 



B. Other pairs of graphs 



Strongly-regular graphs are only one of many classes of 
graphs that are considered difficult to distinguish. Here, 
we discuss two regular (rather than strongly-regular) 
pairs of graphs that are known to be difficult to distin- 
guish (and were, in fact, constructed to be so). The first 
pair, which we denote here as (G\,G2) was first consid- 
ered by Emms et al. [13, [28| (the reader is referred Ap- 
pendix [B] for the adjacency matrices of the two graphs). 
The two graphs are regular having 14 vertices. They were 
given as an example of non isomorphic graphs that can 
not be distinguished by certain quantum random walk- 
ers, although it should be noted that other more recent 
quantum random walkers seem to be able to distinguish 
between the graphs [401 ] . 

While having almost identical adjacency matrices and 
therefore also almost identical energies and Q2 values for 
s — > 1, the two graphs can still be distinguished by our 
quantum annealer, at least within a small region of s. 
This can be seen by looking at differences in Q2 values 
(and smaller but visible differences in the energy). These 
are shown in Fig. [5] Interestingly, this region precisely 
corresponds to the quantum phase transition normally 
seen in quantum adiabatic com put ing procedures that 
is characterized by a small gap [4l|, |42[ ■ This region is 
therefore the usual 'bottleneck' of the Quantum Adia- 
batic Algorithm. 

We also considered another pair of non-isomorphic regu- 
lar graphs on 16 vertices, denoted here by G3 and G4 (the 
reader is referred Appendix [Bl for the adjacency matrices 
of the two graphs) , that are known to be difficult to dis- 
tinguish [3a,|43. Similarly to the previous pair (Gi,G2), 
the pair {G^fii) can be distinguished by looking at the 
Q2 and energy differences as a function of s, see Fig. [6] 

Interestingly, if one inspects the spectra of the 'diagonal' 
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FIG. 6: The difference in Qi values, AQ2 (main panel) and 
energies (inset), for the two graphs G3 and G4 with iV = 16 
that certain quantum walkers cannot distinguish [i^ ]. 

Hamiltonians H p (G) of the two pairs of graphs, one finds 
that each pair shares the exact same spectrum, mean- 
ing that while 'classically' it is impossible to distinguish 
between the graphs within each pair, extending the anal- 
ysis to the 'quantum region', by adding a non-commuting 
term (the driver Hamiltonian) enables the crucial missing 
differentiation. 



V. CONCLUSIONS AND FUTURE RESEARCH 



The results we presented here support a conjecture that 
the Quantum Adiabatic prescription can differentiate be- 
tween all non-isomorphic graphs, given an appropriate 
choice of problem and driver Hamiltonians. This conjec- 
ture needs to be tested more thoroughly, both theoreti- 
cally and also by experiments on real quantum annealers. 

In future studies, it would be desirable to study larger ex- 
amples of SRG's, and also consider other types of graphs 
such as the Cai, Fiirer and Immerman graph construc- 
tions (HI, which are too large to be studied with the 
methods used here and would require Quantum Monte 
Carlo simulations. 

In the current study we have distinguished between 
graphs by measurements of average quantities like M x = 
(^2i a i) (i- e -i magnetization along the ^-direction) and 
J2i^j( a i °j) 2 - More information would be obtained by 



more sophisticated analyses of measurement data. One 
example, would be to calculate higher moments of the ob- 
tained data rather than using only averages as we have 
done in this study. Another example would be compar- 
ing the individual (sorted) values rather than the average 
over sites [45[. This should be done in future studies. 

It is important to understand the efficiency of the algo- 
rithm we propose, i.e. how the running time of the algo- 
rithm scales with graph size for large sizes. This requires 
an analysis of the size dependence of the minimum gap 
for large sizes, presumably using Quantum Monte Carlo 
simulations, e.g. j4f|, to get to large enough sizes to see 
the trend. 

An attractive feature of the method proposed here is that 
it can be implemented in principle by existing (albeit 
prototypical) quantum annealers. We require a quantum 
annealer that has Ising-spin interactions of finite con- 
nectivity plus constant transverse magnetic fields. The 
prototypical "D-Wave One" machine Q has these ca- 
pabilities. However, D-Wave One can not calculate the 
quantities studied here, (of) and (erfcrj), though it can 
calculate averages of a single a z and also susceptibili- 



ties involving a z [47j. For the model studied here, (of) 
is zero because of the bit-flip symmetry of the Hamilto- 
nian. However, one could break this symmetry by adding 
a field coupling to a z and use the values of the (of) to dis- 
tinguish between graphs. In other words, it seems likely 
to us that quantities which D-Wave One can calculate 
would also be useful for the Graph Isomorphism prob- 
lem. Another useful quantity which can be determined 
by D-Wave hardware |47( is the distribution of the energy 
as a function of the evolution time of the algorithm. This 
distribution would be sensitive to the details of spectrum 
of the graph and is therefore also likely be a good dis- 
criminant between graphs. 
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Appendix A: Conjugate gradient method 



In this appendix we explain the version of the conjugate- 
gradient method that was used in this study to calculate 
the ground state of a given Hamiltonian. We found that 
the 'traditional' Lanczos and conjugate-gradient methods 
are not sufficiently accurate for our problem, so we had 
to develop a new variant of these methods. 

The ground state of a Hamiltonian H is obtained by min- 
imizing 



{#>} (w) 



(Al) 



The key feature is that the objective function / has no 
local minima but only the one global minimum that we 
need. Hence, in principle, any minimization method 
would work, but we find that this is not the case in prac- 
tice because of the limitations of finite-precision arith- 
metic. After choosing a basis {!".)}, one can write 
\tjj) = ^2c n \n), and the objective function to be mini- 
mized becomes 



En 



min/({c„}) 

X/ra C ndn + Xm ° n X/r 



(A2) 



mm 



X 



c 2 



in which we denote the diagonal elements of the Hamil- 
tonian by H nn = d n and the off-diagonal elements (pre- 
sumably sparse) by H nm = a nm . In what follows we shall 
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assume that the coefficients c„ are real-valued for conve- 
nience. The generalization to complex- valued coefficients 
is trivial. 

The basic conjugate gradient method in this case is very 
simple. Firstly, the gradient with respect to the various 
parameters c n is easily obtainable: 



df_ 

dc k 



2d k c k + 2J2 m &'krnCm 
V c 2 

2Cfc (Era C ndn + En C n Em a nm c ni) 

(E„4) 2 



(A3) 



Secondly, the one-dimensional minimization steps of the 
conjugate gradient method can be calculated explicitly 



r Q f(\i>)+a\S)) = 

(6\H \S)a 2 + 2 Re{i)\H\5)a + (il>\H\tp) 
1q (6\5)a 2 + 2Rc(ij\S)a + J^tf) 



(A4) 



The above expression is minimized for a = a* which is 
one of the two solutions of a quadratic equation, and this 
produces the minimum (along the line) energy of: 



E a = f(\1>)+a*\S)) 



(A5) 



(5\H\5)a* z + 2Rc{ji\H\5)a* - 

(S\8)a* 2 + 2 Rc(i/j\S)a* + (V#) 



the gradients. Its only effect is to offset the minimum 
energy to zero at each step. The resulting ground-state 
energy at the end of the process is stored in the variable 
offset, Eq = e. 



2. Minimizing the residuals 



Once the energy E is measured to great precision, we 
followed this up with a second stage of minimizing the 
sum of squares of the residuals, where the residuals are 
the elements of H\ip) — Eo\ip). We therefore perform a 
conjugate gradient routine on the objective function: 



._ • (i>\H 2 m~2E (i>\H\^)+E 2 (m 



at the end of which is a good approximation to the 
ground state wavefunction. 



Appendix B: Non-SRG graphs 



Variable-offset minimization 



Because of the limitation of floating point arithmetic, in- 
accuracies may occur when adding terms with very dif- 
ferent orders of magnitude. This turns out to not affect 
the value of ground state energy very much, but it does 
have a large effect on the ground-state coefficients, and 
hence also on ground-state expectation values. To pre- 
vent this, we find it is necessary to ensure that the di- 
agonal and the off-diagonal contributions to the energy 
have the same order of magnitude. 

A simple way to do this is to offset the diagonal term 
in each step such that the total energy is kept equal to 
zero, meaning that the diagonal and off-diagonal terms 
are equal in magnitude (and opposite in sign). This is 
simply done by shifting the diagonal elements d n — > d n + 
e, choosing the constant e in each step in such a way 
that the total energy (diagonal plus non-diagonal) is zero. 
From the equality: 



En C n (d n + <0 + En c ™ E 

EM 2 

En c ndn + En C » Em a nm c n 

E„kl 2 



(A6) 



it follows that offsetting the objective function by e will 
not change the amplitudes c n , and so will also not affect 



Here we provide the adjacency matrices for the models 
studied in Sec. HVB] Both d and G 2 arc regular graphs 
on 14 vertices with valency 4 and are found in [271 128| . 
The graph G\ is given by the adjacency matrix: 



Ac, = 



/0 1 

1 



1 







1 
1 1 
1 1 




1 1 





10 
10 
10 
110 1 
10 
10 1 
10 
10 
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1 
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11 
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10 


1 
10 
110 
1 



\ 1 1 1 1 



1 1 \ 
10 
1 

1 



1 1 
1 
1 

1 

1 

1 

o/ 



(Bl) 



and the graph Gi is obtained from G\ by replacing the 
entries (1,2) and (3,4) with entries (1,3) and (2,4) (and 
corresponding transposed entries) . It can be verified that 
G\ is not isomorphic to G 2 . 



Graphs G3 and G4, also studied in Scc. lIVBl are regular 
graphs on 16 vertices with valency 3 and can be found in 
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Ref. (43[. The graph G3 is given by the adjacency matrix: 
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and the adjacency matrix of G4 is obtained from that 
of G3 by inverting the entries (i.e. interchanging ones 
with zeros) belonging to the sub-matrix spanned by rows 
(1,2) and columns (3,4,5,6) along with the correspond- 
ing transposed entries. 



